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We present an estimate of the nuclear electric polarizability of the 6 He halo nucleus based on 
six-body microscopic calculations. Wave functions are obtained from semi-realistic two-body inter- 
actions using the hyperspherical harmonics expansion method. The polarizability is calculated as a 
sum rule of the dipole response function using the Lanczos algorithm and also by integrating the 
photo-absorption cross section calculated via the Lorentz integral transform method. We obtain 
cue = 1.00(13) fm 3 , which is much smaller than the published value a E ip = 1.99(40) fm 3 [lj ex- 
tracted from experimental data. This points towards a potential disagreement between microscopic 
theories and experimental observations. 
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I. INTRODUCTION 

The nuclear electric polarizability cue is related to the 
response of a nucleus to an externally applied electric 
field. It is an interesting observable because it encap- 
sulates information about the excitation spectrum of a 
nucleus. Recently, it has attracted a lot of attention 
both for light nuclei, see e.g. [l[, and for heavy nuclei, see 
e.g. 0. For light systems the nuclear polarizability is rel- 
evant in the extraction of nuclear quantities from atomic 
spectroscopic measurements. The atomic energy levels 
are affected by polarization of the nucleus due to the 
electric field of the surrounding electrons. Such nuclear 
structure correction, which is proportional to Z 3 o.e/clo 
Q where ao is the Bohr radius, needs to be considered in 
the sophisticated quantum electrodynamics calculations 
of the atomic levels that allow the extraction of charge 
radii from isotope shifts measurements of unstable nuclei, 
see Refs. [3, Q and @ for 6 He and 8 He, respectively. An 
even larger effect of the nuclear structure correction com- 
ing from the polarizability is expected in muonic atoms, 
as the muon mass is larger than the electron mass and 
the orbital radius is smaller. This will be relevant to 
the proposed /x 4 He and /x 3 He experiments Q that aim 
at measuring the nuclear charge radius of 4 He and 3 He 
from the Lamb-shift, to be compared to electron scatter- 
ing data. 

The nuclear electric polarizability of Helium isotopes 
is interesting for the several above mentioned reasons. 
It has been already directly measured or extracted from 
experimental data for the 3,4,6 He isotopes [l[. In the 
case of 3 Hc it is worth mentioning that the data from 
clastic scattering on Pb at energies below the Coulomb 
barrier Q are in disagreement with estimates based on 



calculations of the photo-absorption cross section [l[ , the 
latter being about a factor of 2 smaller. It is also worth 
noticing that the theoretical calculations are in agree- 
ment with photo-absorption experiments, and that the 
band spanned by using different Hamiltonians in the cal- 
culations is smaller than the difference between the data 
taken from photo-absorption cross section and from ion 
scattering experiments. The data analysis involved in 
the latter approach is quite delicate, because one has to 
separate effects of the nuclear force from Coulomb effects. 

In Ref. [|[ the polarizabilities of several Hydrogen and 
Helium isotopes were calculated with ab-initio methods. 
Among the Helium isotopes, 3 He and 4 He were dealt 
with, but no prediction for 6 Hc was provided. It is the 
aim of this paper to fill this gap. 

The 6 He is known as a halo nucleus, where a tightl 
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bound He core is surrounded by two neutrons 
It happens to be the lightest of the known halo nu- 
clei and it is a Borromean nucleus, because the two- 
neutron and neutron-core subsystems are unbound, but 
the three-body system is held together. Due to the very 
small separation energy which characterizes halo nuclei, 
one expects the polarizability of 6 He to be much larger 
than that of the tightly bound 4 He isotope. Experi- 
mental data indicate this behavior. In this paper we 
would like to see whether microscopic calculations repro- 
duce the experimental values and lead to a result where 
a E ( 6 Re) > a( 4 He). 

We perform a microscopic study of the nuclear polar- 
izability a.E for 6 He and compare it to 4 He. We limit 
our study to simple scmi-rcalistic two-body forces. For 
that we use the hyperspherical harmonics method with 
an effective interaction, EIHH, to speed up the conver- 
gence [n], [HJ . The polarizability is calculated as a sum 
rule of the dipole response function using the Lanczos al- 
gorithm and also integrating the photo-absorption cross 
section calculated with the Lorentz integral transform 
method [l3l ]. 

This paper is organized as follows. In Section II we 
describe in details the theoretical calculation of the po- 
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larizability. In Section III we present our results and in 
Section IV we make a comparison with experiment. Fi- 
nally, we conclude in Section V. 



II. THEORETICAL ASPECTS 

The nuclear electric polarizability in the unretarded 
dipole approximation is defined by the expression 
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where a is the fine structure constant, D z is the unre- 
tarded dipole operator and E / n are the energies of the 
nuclear ground and excited states |0) and \n), respec- 
tively. This observable is clearly related to the photo- 
absorption cross section and to the dipole response func- 
tion. The photo-absorption cross section cr 7 (u;) of a nu- 
cleus is given by 



a j (ui) = 47r 2 a d R(u>) 
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where R(uS) is the response function. In the unretarded 
dipole approximation 
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where indicates an average on the initial angular mo- 
mentum projections. The dipole operator is given by 



3 /2, where A is the number of nucleons 



and t£ and Zi are the third component of the isospin op- 
erator and the coordinate of the ith particle in the center 
of mass frame, respectively. One can recover the expres- 
sion for in Eq. ([!} by calculating sum rules of the 
photo-nuclear cross section. The various moments of cr 7 
are defined as 



(4) 
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where ui is the photon energy and u>th and ui indicate 
threshold energy and upper integration limit, respec- 
tively. Assuming that cr 7 (w) converges to zero and utiliz- 
ing the closure of the eigenstates of the nuclear Hamilto- 
nian H, one can relate the polarizability to the n = —2 
sum rule, 
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The polarizability cue can be calculated with the Lanczos 
algorithm using a proper pivot. It is useful to rewrite 
Eq. © as 
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Starting from the "pivot" of Eq. ([7]) where the ground 
state |0) is obtained by solving the Schrodinger equation, 
ole can be expressed as a continued fraction containing 
the Lanczos coefficients XVm, 
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where the \<pi) form the Lanczos 
= 0, . . .}. In fact one has 

oe = -2a(0|£)+ J D z |0) 



orthonormal basis 
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In this work we calculate the polarizability in two dif- 
ferent ways. On the one hand we utilize Eq. ©. On the 
other hand we obtain m_2 by integrating our results for 
the total photo-absorption cross section calculated with 
the Lorentz integral transform (LIT) method [llj]. In 
Ref. [l5|, HU we have presented microscopic calculations 
of the 6 He ct 7 with semi-realistic potential models. Here 
we use lager model spaces which are nowadays available. 
The LIT, an integral transform with a Lorentzian kernel, 
is defined as 
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The LIT is also typically calculated using the Lanczos 
technique explained above, (see [I3|)- In fact it can be 
re-expressed as 
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with z = Eq + a R + iai- It is evident that the LIT 
in (fTT|) is also a continued fraction as in Eq. ((9]), where 
Eq is replaced by a complex z = Eq + a R + iaj. Once 
C((t r ,oi) is calculated, one can invert the LIT to 
get R(lu), and thus m_2. The two methods have to agree 
within the numerical uncertainty. However, with the first 
method one avoids the complications introduced by the 
inversion procedure. 

Given the Hamiltonian H, the calculation of in 
both ways is based on the EIHH [l2[ expansion of the 
wave function. This approach is translational invariant, 
being constructed with the Jacobi coordinates. We use 
different semi-realistic potential models for our calcula- 
tions. Following Ref. (l5j . we will use the Minnesota 
(MN) potential pi] 
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FIG. 1: (Color online) 4 He polarizability: calculations with 
the Minnesota potential for two different u values are shown as 
a function of the grand- angular momentum quantum number 
-Kmax- The polarizabilities obtained from realistic two- and 
three-body interactions 0, HH are presented as a light (blue) 
band. Experimental data from [H. |25| are given by the dark 
(green) band. 



where P^ r are spin and space-exchange operators, Vr,t,s 
are parameterized as linear combination of Gaussians of 
the two-body relative distance and it is a parameter. 
This force reproduces the S— wave nuclcon-nucleon phase 
shifts and correctly binds the deuteron. It renormalizes 
effects of the tensor force into its central component. A 
typical value for u in the Minnesota potential is u = 1, as 
we used in [l5| . Here we will explore the variation of this 
parameter by choosing u > 1. The mixing parameter u 
does not affect the dominant and 3 Si waves in the 
NN interaction, but only affects the s = 1, t = 1 chan- 
nels, where the dominant components are the P— waves 
(■"■Pi and 3 Po,i,2)- For u = 1 there arc no P— waves, they 
contribute only for u > 1. Thus, changing u mostly af- 
fects 6 He, without substantially changing 4 He. Because 
in [11 [II] we also used the Malfliet-Tjon (MTI-III) [1J, 
and the Argonne AV4' [2l[ potentials, we will present 
some results with these interactions as well. The Min- 
nesota potential has been recentl y u sed in a microscopic 
cluster model calculation of 6 He [22| and in the Gamow 
shell model approach [23| for 6 He and 8 Hc. 



III. RESULTS AND DISCUSSION 

The main focus of this work is to study the 6 He po- 
larizability. We start, however, the discussion with the 
4 He nucleus. In Fig. [TJ wc show the results of cal- 
culated via the Lanczos coefficients, as in Eq. ([§]). The 
ground state |0) and the Lanczos "pivot" \<fio) are given 
in terms of the EIHH expansions. While for the ground 
state the expansion is characterized by an even hypcr- 
sphcrical grand-angular quantum number -K m ax and to- 



tal isospin T = 0, T z = 0, D z \0) has to be expanded 
on odd grand- angular quantum number -K max , where the 
isospin in the final state is T' = 1. Fig. Q] shows the 
convergence of ole as a function of -Kmax, where for each 
point -Kmax + 1 is used for the Lanczos "pivot" . We show 
our results for the Minnesota potential with u =1 and 
1.20. The convergence is very good, the dependence on 
u is mild and the results are very close to calculations 
where realistic nuclcon-nucleon (NN) and three-nucleon 
(3N) forces have been used. For the latter one, results 
for effective field theory potentials were presented in @ 
leading to = 0.0683(8) (14) fm 3 (corresponding to the 
upper light (blue) band in Fig. [TJ. The error bar of 
this calculation is accounting for the convergence error 
0.0008 fm 3 and also for the uncertainty in the underly- 
ing dynamics 0.0014 fm 3 . We also show the results of 
ctE for the Argonne vi% two-body force and Urbana IX 
three-body force of Ref. [H[ leading to as = 0.0655(4) 
(corresponding to the lower light (blue) band in Fig. [T]) , 
where the error bar comes from convergence only. The 
experimental data are shown as a darker (green) band. 
These include the more recent evaluation of Ref. [l| based 
on the Arkatov et al. [2(| experimental measurement of 
the photo-absorption cross section and an older result 
reported in Ref. [25II. based on earlier measurements by 
Arkatov et al. [231 . We would like to point out that the 
semi- realistic Minnesota potentials lead to a value of the 
polarizability which is consistent with realistic calcula- 
tions and is only about 15% smaller than the average 
value in the experimental band. 

We can also calculate the polarizability by integrating 
the photoabsorption cross section obtained with the LIT 
method. We get perfect agreement as with the Lanczos 
coefficients. For example, for the standard Minnesota 
potential where u = 1 and for a -K m ax = 12/13 model 
space, the Lanczos method gives = 0.06360 fm 3 and 
integrating a 1 up to 120 MeV we get 0.06336 fm 3 , with 
just a 0.4% difference. 

Now we move to the 6 He nucleus. Wc first calculate 
cxe from the Lanczos coefficients. Also in this case the 
ground state is expanded on even hypersphcrical grand- 
angular quantum number -K m ax, but the total isospin 
is T = 1, T z = —1, and D z \0) is expanded on odd 
^-max = -fr'max + 1 . In this case though, the final isospin 
can be T' — 1 or T" = 2. This leads to two possible 
isospin channels that arc calculated separately and that 
open up at different energies. Experimentally, the T = 1 
channel opens up at photon energy u>th — 0.975 MeV, 
while the T = 2 channel opens up at ojth = 22.77 MeV, 
with 7 6 He — >• 3 Hnnp. Due to the inverse energy weight 
in Eq. ([5]) , the T = 2 channel is expected to be less rele- 
vant to cxe- From our calculations we find that the per- 
centage contribution of the T=2 isospin channel to the 
total polarizability changes from 2% to 4% when varying 
u from 1 to 1.20 in the Minnesota potential. 

In Fig. [21 we present a similar plot as in Fig. Q] for 
6 He with semi-realistic interactions. We observe a much 
slower convergence of cxe for 6 He than for 4 He with all 
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FIG. 2: (Color online) The polarizability of 6 He as a func- 
tion of the grand-angular momentum if ma x for different semi- 
realistic interactions: the Minnesota potential with u = 1 — 1.2 
and the MTI-III potential. 



FIG. 3: (Color online) He two-neutron separation energy 
as a function of the grand-angular momentum K max for the 
Minnesota potential and different u parameters. The experi- 
mental value is also shown. 



TABLE I: Results of the EIHH calculation with K max = 12 
for the different u value of the Minnesota potential. The val- 
ues for the energies are in MeV. 



Potential 


Eq (*He) £ (°He) &n(°He) 


Minnesota 






u=1.00 


-29.949 -30.45 


0.50 


u=1.05 


-29.978 -31.13 


1.15 


u=1.10 


-30.007 -31.88 


1.87 


u=1.15 


-30.037 -32.72 


2.68 


u=1.20 


-30.069 -33.65 


3.59 


MTI-III 


-30.760 -32.24 


1.48 



the potentials employed. By looking at the different u 
values in the Minnesota potential, we sec that the conver- 
gence rate and the value of ole substantially change with 
u. This is related to the variation of the binding energy, 
and consequently of the two-neutron separation energy, 
whose numerical values are shown in Table U for com- 
pleteness. By increasing u we are adding more P— wave 
interactions, which bring additional binding to the 6 He 
nucleus, while leaving 4 He almost unaffected. Naively, 
this makes 6 He more tightly bound and thus more diffi- 
cult to polarize, i.e., ole gets smaller. For the value of 
u = 1 the convergence of the polarizability is particu- 
larly slow, due to the fact that 6 He is barely bound, with 
Sin = 0.56 MeV which is about a factor of 2 smaller than 
the experimental value. With the MTI-III potential we 
get a convergence pattern which is close to the Minnesota 
potential for u = 1.05, because the prediction of Si n is 
similar with these two potentials, see Table HI The infor- 
mation that one gains from Fig. [5] is that by increasing 
Sin we can change the overall slope of the convergence 
pattern of ole- 

For any considered value of u though, it is clear that 



our calculations reproduce the fact that the polarizability 
of the halo nucleus of 6 He is much larger than that of the 
tightly bound 4 He, the ratio being almost an order of 
magnitude. 

Here we would like to point out that Brida and Nunes 
in [22j have used the Minnesota potential with u = 1.15 
in a microscopic cluster model and obtained a separation 
energy Sin = 0.90(5) MeV. This result is different from 
the value we obtain and report in Table |H Their calcu- 
lation is performed without the Coulomb force, but its 
effect cancels in the separation energy. Because for 4 He 
the value reported in [22| for the binding energy is -30.85 
MeV, which is in agreement with our value of -30.86(1) 
(with no Coulomb force), we think that the difference 
is due to the cluster assumption made for 6 He. We do 
not make such assumption and in convergence, the EIHH 
result is exact. In Fig. [31 we show that the separation en- 
ergy Si n is very well converged within the model space 
available for all these potentials. 

We can also calculate the polarizability by integrating 
the photoabsorption cross section obtained with the LIT 
method and compare it to the above results. We quote 
numbers for the Minnesota potential with u = 1.05 in 
the largest available model space K max = 12/13. The 
Lanczos method gives ole =■ 0.7542 fm 3 and integrating 
cr 7 up to 40 MeV (60 MeV) we get 0.7711 fm 3 (0.7827) 
fm 3 . Integrating the cross section we have a 3-4% differ- 
ence, which is due to the fact that LIT is not completely 
converged and the inversion procedure introduces some 
numerical error. 

Now we would like to investigate the dependence of 
the polarizability on the two-neutron separation energy. 
This can be achieved for example by plotting ole versus 
Sin for the different values of the parameter u in the 
Minnesota potential. In Fig. 01 we can see that we find a 
correlation between ole and Si n - Calculations have been 
performed with A' max = 12 (K' max = 13). As an esti- 
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FIG. 4: (Color online) Correlation between oe and £2™ in 6 He 
obtained with the Minnesota potential varying the parameter 
u. The MTI-III and AV4' results are also shown. 



mate of the theoretical error bar in the few-body method 
we take the difference between the largest possible cal- 
culation with if max = 12 and the K max = 8 result. We 
also present the data for the MTI-III and AV4' poten- 
tial (as used in [lfl]) for completeness. The error bars for 
the polarizability increase as the separation energy gets 
smaller. This is a reflection of the slower convergence 
observed in Fig. [2j For the Minnesota potential S^n has 
a negligible error, hardly visible in Fig. 21 For the MTI- 
III and AV4' potentials the error in S^n is large because 
these interaction models are not as soft as the Minnesota 
force. 

In Ref. [j| it was argued that the polarizability should 
roughly scale like the inverse square of the binding energy 
of a nucleus. For a halo system, like 6 He, the relevant 
scale parameter is the separation energy, rather than the 
binding energy The - S^n dependence empirically 
observed in Fig- [4] is compatible with such a behavior. 

In order to reproduce the polarizability of a halo nu- 
cleus, it is expected that the halo structure, thus S^n, 
should be correctly modeled, even if the absolute binding 
of 4 He and 6 He are not reproduced. Thus, one can esti- 
mate the value of «£ by choosing S^n to be around the 
experimental value and then calculate the corresponding 
polarizability. A value of u that gives S^n close to exper- 
iment is u = 1.05, where the convergence of cue is slower 
than for larger values of u. From a closer look to Fig. [5] 
and Fig. [3] we can see that also for u = 1.20 the polariz- 
ability oiE is still increasing when K max becomes larger, 
even though the separation energy is converged. This 
means that the convergence of the polarizability is not 
only influenced by S^n- Other observable that is natu- 
rally related to the polarizability in the unretarded dipole 
approximation is the radius operator. 

In recent papers [HI, [2^] the correlation between the 
polarizability and the neutron skin of the 208 Pb nucleus 
was studied within the nuclear density functional theory 
framework. In the following, we will investigate the same 
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FIG. 5: (Color online) Correlation between the nuclear elec- 
tric polarizability of 6 He and the skin radius with the Min- 
nesota potential with different u. The three points for each 
u value correspond to calculations with K max = 8, 10 and 12, 
moving from the left to the right. 
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FIG. 6: Neutron skin radius r s ki n of 6 He with the Minnesota 
potential and it = 1.05, as a function of the grand-angular 
momentum quantum number if max . The curve is a fit to the 
calculated points, used to extrapolate to infinite model space. 



correlation for 6 He, even though 6 Hc is a different system. 
For halo nuclei, one refers to the halo radius, rather than 
the skin radius, but clearly the observable 



^*skin 



(13) 



where r n and r p arc the mean point-neutron and point- 
proton radii, can be uniquely defined. In our recent 
work J^jj, such observables have been calculated for 6 He 
from realistic two-body potentials in the EIHH method. 
Here, instead, we use the same semi-realistic interaction 
as for the cxe calculations. In Fig. [5j we show a plot of 
cue versus r s ki n for different model spaces and for three 
different values of u in the Minnesota potential. The 
three points for each u correspond to calculations with 
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K max = 8, 10 and 12, from the lowest to the largest value 
of Oe, respectively. We clearly see a linear dependence 
between a e and r s ki n for all the three u values where the 
coefficients depend on the separation energy as 



a E = a(S 2n ) + b(S 2n )r t 



kin - 



(14) 



Because S2 n is converged and because of the linear de- 
pendence displayed in Fig. [5] we deduce that the calcu- 
lation of ole is not fully converged because the radii, 
and especially r„, are not fully converged. The calcu- 
lation of a radius of the ground state does not require 
an expansion on the dipole excited states as in Eq. ([7)1 
and as such is less computation demanding and can be 
performed for larger model spaces, where radii arc bet- 
ter converged. Thus, the approach we take to estimate 
a E from our calculations is to "measure" the coefficients 
a(S2n) and b(S2n) from the ci£ results in the available 
model spaces and, assuming that this physical linear de- 
pendence will be unchanged in larger model spaces, we 
will use the coefficients to obtain aE from a bound state 
calculation of just r sk in- Choosing S2n close to exper- 
iment, in the case of u = 1.05, we find a = —1.275 
fm 3 and b = 2.616 fm 2 . We then calculate r sk j n , see 
Fig. |6l up to the largest grand-angular momentum value 
accessible with our computational facility, A max = 16. 
One can clearly see that convergence is approached. Ex- 
trapolating these points with an exponential ansatz of 
the form r skin (A" max ) = r skin (oo) - ce^^- we get 
r s kin(oo) = 0.87(5) fm. As an error estimate we take the 
difference between r skin (A max = 16) and r skin (A max = 
12). The theoretical value is somehow larger than the 
experimental data. In fact, combi ning different mea- 
surements of the matter radius [3ll433[ with the most 
recent evaluation of the proton radius [5J, one can infer 
r n and consequently the skin radius, which is found to 
be 0.52< r^P <0.62 fm. The variation on r^f n is fairly 
large, due to the large uncertainty in the matter radius 
determination from ion scattering. 

At this point, from our extrapolated skin radius, using 
Eq. (fl4"|) we estimate the theoretical nuclear electric po- 
laziability of 6 He to be = 1.00(13) fm 3 . If we were 
to use Eq. (|T4"f with the experimental values of the skin 
radius, one would obtain a nuclear electric polarizability 
of aE = 0.08 — 0.34 fm 3 , which is even smaller than the 
estimate based solely on theory. 



IV. COMPARISON WITH EXPERIMENT 



In Ref. [l| an evaluation of the experimental number 
for the polarizability of 6 He was presented, leading to 
1.99(40) fm 3 . This was obtained from the inverse 



cxp 



energy weighted integral of experimental and theoreti- 
cal B(E1) response functions for 6 Hc. The experimental 
distribution was measured by Aumann et al. [34j from 
the Coulomb breakup of 6 He on Lead and Carbon at 




5 10 15 20 25 30 35 40 45 50 
co/co [MeV] 

FIG. 7: (Color online) Double inverse energy weighted 
cross section as a function of the energy u and sum rule 
m-2(oj)/2ir 2 as a function of u) for 6 He with the Minnesota 
potential and u — 1.05. 
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FIG. 8: (Color online) Double inverse energy weighted cross 
section as a function of the energy u for 4 He and 6 He with 
the Minnesota potential and u = 1.05. 



240 MeV/u up to 8 MeV above threshold. In order to 
obtain the polarizability, data were extrapolated up to 
12.3 MeV, where the threshold for the breakup into two 
tritons opens up. The theoretical curve was taken from 
a calculation of the dipole transition to the 1 _ contin- 
uum (3oT | in a three-body model with phenomenological 
n — n and n — a interactions plus an effective three-body 
force, also calculated up to the two tritons threshold. 
The estimate was done in two steps: (i) an average of 
the experimental data and theoretical curve was taken 
up to 12.3 MeV; (ii) to account for the higher energies, 
the polarizability of 4 He was added. The latter one basi- 
cally comes from integrating the photo-dissociation data 
from Arkatov et al. [26| . Because we can access the full 
response functions using the LIT method at any energy 
below the pion production threshold, we can verify these 
two approximations. First, it is interesting to study the 



7 



convergence of ole calculated as a sum rule of the re- 
sponse, see Eq. (0} and (JSJ) , to investigate the validity 
of (i). In Fig. [3 we present both the integrand func- 
tion cr 7 (w)/27r 2 cj 2 versus uj and the convergence of the 
integral m_2(w)/27r 2 versus Q. At ui = 8 MeV the sum 
rule is exhausted only up to 75%. Thus only 75% of the 
a E ip is based solely on experimental data. At U) =12.3 
MeV, where the two 3 H channel opens up the sum rule 
is exhausted up to 90%. We observe that one needs to 
integrate up to 40 MeV of energy to have the sum rule 
exhausted at the 98% level. 

To verify the approximation (it) we can compare the 
integrand function ct 7 (w) /2tt 2 ui 2 for 6 He and 4 He at ener- 
gies beyond 12.3 MeV. In Fig. flHJ we observe that the two 
curves agree with each other for lu > 35 MeV, where the 
sum rule is almost exhausted. In the region beyond the 
4 He disintegration threshold and below about 35 MeV 
one would over estimate the sum rule integrating the He 
curve, because one gets 0.044 fm 3 , to be compared to the 
0.027 fm 3 obtained when correctly integrating the 6 He 
curve. On the other hand, neglecting the part of the 
cross section for ui > 12.3 MeV and below the 4 He dis- 
integration threshold one underestimates the sum rule. 
The contribution of this portion is 0.037 fm 3 , about 5% 
of the sum rule. These two effects almost cancel out so 
that, the approximation (it) does not lead to a big error. 

We think that the main reason of the disagreement 
between the estimate from Ref. [l| and our calculations 
comes from the difference in the low-energy part of the 
response. In our previous work [l5l IT(| we have shown 
that our calculations with semi-realistic potentials un- 
derestimates the data from Aumann et al. 34j. Thus, 
what we observe for the polarizability is consistent with 
this fact. We would like to point out that (i) nuclear 
corrections mi ght affect the results in the ion scattering 
experiment of [34[ and that (ii) as discussed earlier, sim- 
ilar experiments for 3 He lead to large discrepancy with 
photo-dissociation results. Nevertheless, measuring ole 
from the dipole response function it would be desirable 
to have data up to higher energies than 12.3 MeV. Addi- 
tional or alternative measurements of oie would help to 



better contraint this observable. 



V. CONCLUSIONS 

We summarize our results as follows. We have carried 
out an estimate of the nuclear polarizability of 6 Hc based 
on the hyperspherical harmonics expansion with simple 
semi-realistic potentials. Our calculations clearly repro- 
duce the fact that the polarizability of the halo nucleus of 
6 He is much larger than that of the tightly bound 4 He. 
For 4 He the semi-realistic Minnesota potentials lead to 
a value of the polarizability which is consistent with re- 
alistic calculations and is about 15% smaller than the 
average value in the experimental band. Nevertheless, a 
large disagreement is found for 6 He. In order to estimate 
cue we have chosen a potential that reproduces the sepa- 
ration energy and then we investigated the correlation of 
the polarizability with the skin radius. Our final result is 
ole = 1.00(13) fm 3 , which is about a factor of 2 smaller 
than the estimates from experimental data. This points 
towards a disagreement of microscopic theory and exper- 
iments. To shed light on this, it would be nice to have 
more data or alternative measurements of ole- Concern- 
ing the theoretical calculations, it is desirable to extend 
these results to realistic potentials including also three- 
body forces. We leave this subject to a future work. 
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